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It was shown recently that a PT-symmetric iip quantum field theory in 6 — e 



. dimensions possesses a nontrivial fixed point. The critical behavior of this theory 

\ around the fixed point is examined and it is shown that the corresponding phase 

transition is related to the existence of a nontrivial solution of the gap equation. 
' The theory is studied first in the mean-field approximation and the critical exponents 

are calculated. Then, it is examined beyond the mean-field approximation by using 
renormalization- group techniques, and the critical exponents for 6 — e dimensions are 
calculated to order e. It is shown that because of its stability the PT-symmetric icj)^ 
theory has a higher predictive power than the conventional (j)^ theory. A comparison 
' of the icf)^ model with the Lee- Yang model is given. 
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^ ■ I. INTRODUCTION 

^ \ The study of 'PT-symmetric quantum theory [1| was originally motivated by the discovery 

that the eigenvalues of the quantum-mechanical Hamiltonian H = p^ + ix^ are real, positive, 
and discrete j2[. This work naturally led to a number of studies of the properties of a scalar 
quantum field theory having an imaginary cubic self-interaction term |3[. A very recent 
study of the renormalization group (RG) equations for the PT-symmetric igcj)^ quantum 
field theory in c? = 6 — e dimensions shows the existence of a nontrivial fixed point . This 
allows for a nonperturbative renormalization of the theory and suggests that the theory 
undergoes a continuous phase transition. 

In the present paper we study this transition in detail; that is, we examine the critical 
behavior of the theory. We will see that such a transition is associated with the existence 
of a nontrivial solution to the gap equation at a critical value ml of the bare mass m?. The 
correlation length (the inverse of the renormalized scalar mass) diverges at im? = ml . 

In Sec. HI] we study this transition within the framework of the mean- field approximation 
and in Sec. Illll we obtain the critical behavior of the theory by calculating the critical expo- 
nents. Next, in Sec. IIVI we compare the conventional and the PT-symmetric 0^ theories and 
analyze the relation between their renormalization and stability properties. We show that 
compared with the conventional 0^ model, the PT-symmetric icfr" theory exhibits new and 
interesting features. Remarkably, the PT-symmetric theory has a higher predictive power 
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than the (f)^ theory. This is because its critical behavior (and therefore its renormahzation 
properties) are governed by one parameter less than the conventional theory. We will show 
that this property is related to the different stability properties of the two theories. 

Section |V] goes beyond mean-field analysis. We study the theory near = 6 dimensions 
and calculate the critical exponents up to 0(e). For (i < 6, the fluctuations around the 
mean-field configuration become important and the analysis of the critical behavior requires 
the use of RG techniques. With the help of the hyperscaling relations, the critical exponents 
are calculated for (i = 6 — e. Some conclusions are given in Sec. |Vll 

II. MEAN-FIELD ANALYSIS 

In Ref. j3] a nontrivial fixed point of the RG equations of the "PT-symmetric igcj)^ quantum 
field theory in d = 6 — e dimensions was found. The existence of such a fixed point suggested 
the onset of a continuous transition. In this section we study this transition by performing 
a mean-field analysis. 

The partition function Z[h] of the quantum field theory in d dimensions is given by 

= y"p(/.e-^[*l"'^'''"^^ (1) 

where h is an external field. The action is S[(f)\ = J d'^x [(9^0)^/2 + V^(0)] and the potential 
is V{(t)) = m2072 + ig(f)^/6. 

We consider the mean-field approximation to Z[h], first searching for a constant-field 
saddle point and then performing a semiclassical expansion around this configuration. [The 
exponential in ([T]) contains an implicit factor oil/h, and in the semiclassical approximation 
h is treated as small, which justifies the use of steepest-descent asymptotic techniques [5|.] 
Without loss of generality, we assume that g > 0. 

The constant-field saddle points are given by the gap equation 

m'^(f) + ig(j)'^/2 = -ih. (2) 

When h = 0, becomes 

0(m2 + ig^/2) = 0, (3) 

which has the two solutions = and 02 = 2im? / g. We must determine which of the 
saddle points, or 02? contributes to the small- /i asymptotic behavior of the path integral 

A. Asymptotic analysis for d = 

Let us first examine the problem for the simple case d = 0. The second derivative of 
the exponential in ([T]) at /i = is —m? — igcf). Thus, at = 0i the second derivative is 
— m^, and at = 02 the second derivative is m^. The sign of the second derivative of the 
potential determines the directions of the steepest-descent paths (constant-phase contours) 
in the neighborhood of the saddle points. Thus, if m? is positive, the down directions from 
01 are horizontal (parallel to the Re0 axis) and the down directions from 02 are vertical 
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(parallel to the Im0 axis). However, if is negative, the directions are reversed; the down 
directions from 02 are horizontal and the down directions from (pi are vertical. 

It is necessary to deform the original integration path in ([1]), which lies on the real 
axis, into a constant-phase contour that passes through the appropriate saddle point. If 
we let = M + iv and note that the phase (the imaginary part of the exponent) at the 
saddle points vanishes, we obtain a cubic algebraic equation for the constant-phase contour: 
u {gv'^/2 — w?v — gu^/6) = 0. Thus, near oo the steepest-descent contours asymptote to the 
angles 7r/2, — 7r/6, and — Svr/G and the steepest-ascent paths asymptote to the angles — 7r/2, 
7r/6, and 5n/6. Thus, the uniquely determined steepest-descent contour is a PT-symmetric 
(left-right-symmetric) path that terminates at the angles —ir/G and —bir/G. 

Hence, the crucial observation is this: When > 0, 02 lies above 0i on the Im0 axis 
and the steepest-descent path passes through 0i (and not 02)- However, when < 0, 02 
lies below 0i on the Im axis and the steepest-descent path passes through 02 (and not 0i). 
A phase transition occurs at = where the saddle points coincide. 



B. Asymptotic analysis for d > 

When d > 0, we consider the Hessian matrix D~^{q,p) = A~^(g^)5'^(p + q) in Fourier 
space, where 

= q^ + + ig(f)i. 

This is just the inverse tree-level correlator when 0j is the vacuum. For > 0, D'^{q,p) 
is positive definite when evaluated at 0i = 0. In this case we find that 

A'\q^) = q^ + m\ (4) 

On the other hand, for < 0, D~^{q,p) is positive definite when calculated at 02 = 2iiv? / g, 
and we get 

A-^(g2) = g2 -2m2 = g2 - m^. (5) 

These results imply that there is a phase transition in the Euclidean "PT-symmetric igcj)^ 
theory, with the two phases being determined by the order parameter 0, where 0i = and 
02 = 2im^ / g. 

We proceed to identify the relevant couplings-^ in general, these are the couplings that 
control the phase transition. In the ferromagnetic case, the relevant couplings are the tem- 
perature and the external magnetic field j6|. Thus, we look for the two couplings that in 
our case play the role that the temperature and the external magnetic field play in the 
ferromagnetic case. We begin by defining the temperature T as 

T = 2m' /g, (6) 

where the coefficient 2/ g in (jS]) is chosen for later convenience [3]. The critical temperature 
Tc is obtained with the help of ([3]) for the nontrivial saddle point (that is, the solution to 
m? + ig(t)/2 = 0). We seek the limiting value of T for which vanishes, and in our case we 
have Tc = 0. Therefore, the reduced temperature is given by 



T = T-T, = 2m' /g. 



(7) 
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The external magnetic field (the conjugate field) is just the parameter h. 

From (jl]) and we see that the excitations above the trivial vacuum 0i = and those 
above the nontrivial one 02 = ^irn^/g have the same mass, M = |m^|^/^, which goes to 
zero as we approach the critical point. Moreover, as we move continuously from positive to 
negative values of m^, = 2iim? / g moves continuously down the negative-imaginary axis 
starting from = (when m? > 0). The presence of such a divergent correlation length 
(the inverse of the scalar mass M = 0) and of a continuously varying order parameter are 
the indications of a continuous (second-order) phase transition. In RG language, this is due 
to the existence of the nontrivial fixed point found in Ref. 0]. 

We have shown that at least within the framework of the mean-field approximation used 
in this section, the fixed point found in Ref. {3] governs the transition from the 0i = phase 
to the 02 = 2im? / g phase. In Sec. IIIII we calculate the critical exponents that govern the 
behavior of the theory in the critical region in terms of the parameters h and r. 



III. CRITICAL EXPONENTS IN THE MEAN-FIELD APPROXIMATION 

In this section we give a quantitative description of the critical behavior of the model 
for d > Q dimensions. We define and calculate the critical exponents in the mean-field 
approximation by following the standard terminology used for Ising-like models. With the 
help of ([7]) we write the saddle-point equation (|2]) for /i 7^ as 

r + z0 = -2zV (8) 

which has the typical form of a gap equation. From (|8]) we can immediately deduce the 
mean-field exponents /3 and 5. The exponent /3 is determined by 0|^_q ~ t^, which gives 
(5 = 1. The quantity 5 is defined by the relation 0|^_q ~ h^^^ , which implies that 5 = 2. 

To evaluate the exponents 77, and 7, we need the two-point inverse correlator given 
in (j4]) and (|5]) for the = and the = 2im? / g phases, respectively. The exponent rj is 
defined by A(g^) ~ q^""^ (g^ — t- 00), which implies that 77 = in both phases. 

To calculate v we need the correlation length. The latter is proportional to the inverse 
of the renormalized mass M (the pole of the propagator). In the mean-field approximation 
discussed in the present section, = \m?\. For the phase with nonzero (0 = 02) we find 
that = —m? ~ 02 ~ r. This implies that the mean-field exponent which is defined 
by ^ ~ r-^ is 1/ = 1/2. 

Finally, the exponent 7 is obtained as follows. The susceptivity x is given by x = 
i 50/5/i|^^Q. By differentiating ([2]) with respect to h we have m'^6(j)/6h + ig(j)5(j)/5h = —i, 
which we evaluate at /z = to get (m^ -|- igcp) x = 1- We thus obtain 

X = (m^)"! ~ for T > T, (m^ > 0), 

X = {-m^)-^ ~ for T < Te (m^ < 0). 

Therefore, 7, which is defined by the equation x ~ t ''^, turns out to be 7 = 1. 

In this and in the previous section we have studied the features of the continuous tran- 
sition of the PT-symmetric ig(j)^ theory in the mean-field approximation by calculating the 
critical exponents. In the following section we deepen our understanding of the behavior 
of the theory in the critical region. To this end, we compare the conventional and the 
PT-symmetric theories by relating their critical behaviors to their stability properties. 
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IV. RENORMALIZATION AND STABILITY 

In Ref. [3] we compared the conventional (f)^ theory with the corresponding icf)^ VT- 
symmetric theory near and at ci = 6 dimensions, and the instabihty of the former was 
contrasted with the stabihty of the latter. Moreover, for d = Q — e dimensions, the perturba- 
tive renormalizability of the 0^ theory (obtained by taking the bare parameters around the 
Gaussian fixed point) was compared with the renormalizability of the PT-symmetric theory. 
The latter is nonperturbative because it is realized around a non-Gaussian fixed point (which 
collapses onto the Gaussian fixed point at d = 6). 

As shown in Ref. [3], there is a crucial difference between the RG properties of the conven- 
tional and the PT-symmetric theories. In the former both and g are relevant directions, 
but in the latter is the only relevant direction {g being irrelevant). Accordingly, in 
Sec. mil we found that to make the nontrivial solution 02 vanish we only need to tune the 
temperature T = 2m? / g towards its critical value Tc = 0. For generic values of the bare 
coupling g (which is kept fixed), we must tune only im? towards zero. This is similar to what 
happens in conventional 0^ theory, where the only relevant direction is im? and the quartic 
coupling plays no role in determining the phase transition. It appears that the instability 
of the 0^ theory gives rise to an additional condition to reach the critical region. 

In this section we examine the origin of the profound difference between the two theories. 
To this end, we study the mean-field (classical) potential of the conventional theory 

y(0) = mVV2 + ^0V6- (9) 

As in the previous section, without loss of generality we consider the case (7 > 0. 

This potential is unstable for both the m? > and < cases (see Fig. [1]). However, if 
the false vacuum is sufficiently long lived (that is, if the tunneling time is sufficiently long), 
the theory can be consistently defined around this vacuum. For h = the gap equation 
for the potential in is (pirn? + g(j)/2) = 0. For > 0, the potential has a minimum 
at the origin = and a maximum at = —2m?/g\ for m? < 0, there is a minimum at 
= —2m?/g and the maximum is at = 0. Defining the reduced temperature r as 

r = w?/g, (10) 

we study the continuous phase transition from = —2m?/g to = 0. We stress that in 
the "PT-symmetric igcfr' theory (as well as in the conventional (70^ theory) the scaling region 
r — is reached by keeping g finite and taking the limit m? — >■ 0. This is not true for the 
ordinary 0^ theory. 

The potential in ([9]) for m? < is given by the solid line in Fig. ([1]) [apart from a trivial 
constant shift ^(0) — )■ ^(0) — 2m^ / (3(7^)]. In the limit rr? — > with g = finite the potential 
becomes the dotted line (cubic parabola), which has an inflection point at = 0. In this 
limit the vacuum disappears. 

Let us define the area A, 

Jrrfi/g °9 

of the surface under the potential function between the points = —2m?/g and = m?/g 
(see Fig. [T]). Note that in the limit — )■ with g = finite the nontrivial minimum 
= —2rn?/g moves toward = 0, so ^ decreases. Therefore, the lifetime of the false 



(11) 



6 



vacuum, which is proportional to A, becomes smaller, and the theory is destabilized. As 
already mentioned, from the RG viewpoint this happens because in addition to m^, g is a. 
relevant parameter. To reach the continuous phase transition (r — 0) while keeping the 
theory stable {A = fixed), we must tune not only but also g. From ( ITO|) and (fTT]) we get 

~ A/t^ and g ~ A/t"^. (12) 

This means that in order to reach the continuum limit, and 1/g must be tuned sepa- 

rately to zero with r — ?■ according to (fT2|) while keeping A (the vacuum lifetime) fixed. It 
is clear that the conventional 0^ theory has an additional relevant direction as compared to 
its PT-symmetric counterpart because of the intrinsic instability of the theory. The VT- 
symmetric model, being energetically stable, has only one relevant direction, thus showing 
a higher degree of predictive power. 
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FIG. 1: The potential V{(j)) = m'^cf? /2 + g<jr' /Q-2m^ / {'ig'^) for g = 6 x lO'^ and = -1.8 x lO'^ 
(solid line), = —1.2 x 10~^ (dashed line) and = (dotted line). For the solid and the dashed 
lines, the vacuum is at = —2m?'/g while the intersection with the negative axis \s cf) = rr? jg. 
The area A of the surface included between the ^ axis and T^(</') is proportional to the lifetime of 
the false vacuum. 



V. CRITICAL EXPONENTS NEAR d = 6 

In this section we study the critical behavior of the theory for c? < 6. (The mean- field 
results provide a good description for > 6.) We turn our attention to the calculation of 
the critical exponents of the "PT-symmetric theory beyond the mean-field approximation 
considered in Sec. Illli [According to the mean-field analysis of Sec. Illlt in the presence of 
the external source h the two relevant parameters are t = T — T^. = 2m? / g (the reduced 
temperature) and h itself (the external field) .] As is known from the general theory of critical 
phenomena, below the upper critical dimension (d = 6 in this case), the fiuctuations around 
the mean-field configuration became important. RG techniques provide an essential tool for 
calculating the scaling behavior of the theory. (Note that although the exponents rj and 
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5 appear in the work of Fisher |8[ on the Yang-Lee zero problem for the first time, the 
evaluation of the other exponents is an original achievement of this paper.) 

Let us consider our theory in^d = 6 — e dimensions. The RG equations imply the existence 
of two nontrivial fixed points 



h* = 0, m^* = 0, g* = ±V1287r3e/3. 

To each of these points is associated the phase transition that we have just described in 
the mean-field approximation. (The above analysis was done for g* > but for g* < the 
results are analogous.) 

According to Ref. [4] the scalings of h, m?, and g with the running scale t = ln(/i/ /iq) are 
given by h{t) = cie^^*, m^(t) = C2e^^\ and g{t) = g* + c^e^^*, where 

^i = _4 + 4e/9, (72 = -2 + 5e/9, ^73 = 6. (13) 

Prom the above equations we see that the two relevant parameters are h and r = 2m? / g. 
With the help of the hyperscaling relations 

o d + gi gi 2gi + d 1 0,^,0 

P = , ^ = ' 7 = , ^ = , r] = 2 + d + 2gi, 

92 d + gi g2 g2 

we can calculate from (fT3|l the critical exponents, which turn out to be 

[3 = 1, (14) 

7 = l + e/3, (15) 

6 = 2 + e/3, (16) 

u = l/2 + 5e/36, (17) 

r] = -e/9. (18) 

From ( !T8|) . the scaling dimension of the scalar field is [0] = (c/ — 2 + ?7)/2 = 2 — 5e/9. Note 
that for e = these exponents coincide with the mean-field values calculated in Sec. IIL 



VI. CONCLUSIONS 

The model was introduced in ^ to study the density of the Lee- Yang zeros of the 
partition function Z[H] on the imaginary axis of H [H being the magnetic field). The 
"critical exponents" in Ref. [sl are not critical exponents in the physical sense, but rather 
they are parameters governing the mathematical behavior of the function that gives the 
asymptotic density of zeros on the H imaginary axis near the branchpoint H = 0. If we set 
rr? = in our theory, we obtain the Lee- Yang model studied in Ref. |8|, where the exponents 
rj and 6 already appear; /3, 7, and u are given for the first time in (fT4|) . ( !T5|) . and (fTTl) . 

In contrast to the Lee- Yang model, our ig(j)^ theory has a physical interpretation as a VT- 
symmetric Euclidean quantum field theory. Moreover, to show that our theory undergoes 
a second-order phase transition at = (that is, to show the renormalizability of the 
theory) we have had to investigate its behavior in the neighborhood of = 0; it was not 
sufficient to study the theory at = 0. We emphasize that the physical symmetry of the 
theory (that is, the VT symmetry) allows for the presence of the 0^ operator in addition to 
the linear and cubic terms. 
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Using RG techniques, we have studied the theory beyond the mean-field approximation 
and have calculated the critical exponents for d = 6 — e dimensions up to 0(e). In studying 
the critical behavior of the PT-symmetric ig(j)^ quantum field theory, we have shown that the 
phase transition is associated with the existence of a nontrivial solution of the gap equation 
at a critical value ml of im? . We conclude that one can view the Lee- Yang model considered 
in Ref. jsl as the critical theory of the "PT-symmetric igcj)^ model. 

Compared with the conventional 0^ model, the PT-symmetric theory exhibits new and 
interesting properties. In particular, it has a higher degree of predictive power because its 
critical behavior is governed by one parameter less than in the 0^ theory. We have shown 
that this crucial difference is related to the different stability properties of the two theories. 
Thus, from our work it appears that the renormalization properties of the PT-symmetric 
igcjr' model, when compared with those of the conventional cfr" theory, are quite remarkable 
and encouraging for further studies and future apphcations. 
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